IDA 


INSTITUTE  FOR  DEFENSE  ANALYSES 


Impulse  Response  of  Alternative 
Synthetic  Apertures  for 
Subsurface  Detection 


James  M.  Ralston 


June  2001 

Approved  for  public  release; 
distribution  unlimited. 


IDA  Paper  P-3598 
Log:  H  01-000409 


Report  Documentation  Page 


Report  Date  Report  Type 

01JUN2001  Final 

Dates  Covered  (from...  to) 

01JAN1999  -  01  OCT  1999 

Title  and  Subtitle 

Contract  Number 

Impulse  Response  of  Alternative  Synthetic  Apertures  for 
Subsurface  Detection 

Grant  Number 

Program  Element  Number 

Author(s) 

Project  Number 

Ralston,  James  M. 

Task  Number 

Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es) 

Institute  for  Defense  Analyses  1801  N.  Beauregard 

Street  Alexandria,  Va  22311-1772 

Performing  Organization  Report  Number 

Sponsoring/Monitoring  Agency  Name(s)  and 
Address(es) 

Sponsor/Monitor’s  Acronym(s) 

Institute  for  Defense  Analyses  1801  N.  Beauregard 

Street  Alexandria,  Va  22311-1772 

Sponsor/Monitor’s  Report  Number(s) 

Distribution/Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

The  original  document  contains  color  images. 

Abstract 

Subject  Terms 

Report  Classification 

unclassified 

Classification  of  this  page 

unclassified 

Classification  of  Abstract 

unclassified 

Limitation  of  Abstract 

SAR 

Number  of  Pages 

22 

This  work  was  conducted  under  IDA’s  independent  research  program, 
Central  Research  Project  CRP-2057.  The  publication  of  this  IDA 
document  does  not  indicate  endorsement  by  the  Department  of  Defense, 
nor  should  the  contents  be  construed  as  reflecting  the  official  position  of 
that  Agency. 

©  2001  Institute  for  Defense  Analyses,  1801  N.  Beauregard  Street, 
Alexandria,  Virginia  22311-1772  •  (703)  845-2000. 

This  material  may  be  reproduced  by  or  for  the  U.S.  Government. 


NSTITUTE  FOR  DEFENSE  ANALYSES 


IDA  Paper  P-3598 


Impulse  Response  of  Alternative 
Synthetic  Apertures  for 
Subsurface  Detection 


James  M.  Ralston 


PREFACE 


This  paper  was  prepared  under  an  IDA  Central  Research  Project  entitled  “Impulse 
Response  of  Synthetic  Aperture  Radar.” 
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I.  INTRODUCTION 


A.  BACKGROUND  AND  MOTIVATION 

The  application  of  radar  to  the  detection  and  identification  of  subsurface  targets  is 
generally  a  struggle  to  distinguish  the  targets  sought  from  manmade  and  naturally 
occurring  clutter  objects  [Refs.  1,  2,  3].  Radar  imaging  [Ref.  4],  which  has  been  under 
development  by  many  organizations,  is  one  method  of  achieving  this.  Several  approaches 
to  radar  imaging  have  received  attention,  including  microwave  holography,  tomography, 
and  synthetic  aperture  radar  (SAR).  The  common  element  to  all  of  these  is  the 
measurement  of  the  scattered  radiation  from  target  and  background  from  several  distinct 
locations  in  real  space,  each  of  which  illuminates  the  scene  from  a  different  aspect  angle. 
The  data  from  these  measurements  is  then  analyzed  by  any  of  several  different  transforms 
or  inverse  methods  to  arrive  at  an  image  of  the  target+background.  The  objective  of  the 
effort  is  to  obtain  a  sufficiently  detailed  image  that  the  targets  sought  can  be  distinguished 
from  clutter. 

B.  ILLUMINATION  FUNCTION  AND  IMPULSE  RESPONSE 

All  of  these  imaging  techniques,  whether  two-  or  three-dimensional,  fundamen¬ 
tally  proceed  by  illuminating  the  target  region  with  a  wideband  waveform  from  a  range  of 
aspect  angles  or  directions.  This  illumination  can  be  characterized  as  an  illumination  or 
transfer  function  defined  over  a  volume  of  k-space  [Ref.  4], 

p(k)=  J drp(r) exp(2 jk  ■  (R  +  r ))  ,  ( 1 ) 

where  k is  the  wavenumber  vector  defining  the  direction  and  frequency  of  a  plane  wave, 
p{  r  )  is  the  scattering  density  of  the  target  region  being  imaged,  R  is  the  vector  between 
the  radar  and  the  object  being  imaged,  andr  spans  the  volume  of  integration.  If  the 
scattering  density  is  a  spatial  impulse  Sir  ),  and  the  vector  R  is  set  equal  to  zero,  which 
can  be  done  without  loss  of  generality,  Eq.  1  can  be  inverted  to  obtain  an  estimate  of  the 
spatial  impulse  response  of  the  system 


where  the  summation  is  taken  over  the  set  of  /c- vectors  used  to  illuminate  the  target 
region,  and  TV  is  a  suitable  normalization  factor.  This  estimate  is  equivalent  to  the  impulse 
response  of  the  system  under  consideration. 


C.  CROSS-RANGE  RESOLUTION  AND  RANGE  OF  ASPECT  ANGLE 

Equation  2  is  sufficiently  general  to  apply  to  all  cases  of  microwave  imaging.  A 
useful  expression  for  the  cross-range  spatial  resolution,  Ax,  of  an  imaging  system  is 

.2/1 

Ax  = - = -  ,  (3) 

4sin(A<9/2)  2-A0 

where  2  is  the  wavelength  of  the  illumination  and  A6  is  the  angle  subtended  by  the 
synthetic  aperture  at  the  imaging  point.  Equation  3  states  that  the  resolution  in  the  direc¬ 
tion  perpendicular  to  an  illumination  direction  is  inversely  proportional  to  the  apparent 
angular  extent  of  the  synthetic  aperture  normal  to  the  illumination  axis.  When  the 
imaging  point  is  below  the  surface,  however,  the  angle  is  compressed  by  surface 
refraction  in  the  following  way: 


By  Snell’s  law: 


sin(<9;. )  =  n2  sin(<9r )  ; 


(4) 


for  angle  differentials: 

dOr  _  «i  cos(6>.) 
d6i  n2  cos(0r ) 

where  n\  and  ni  are  the  refractive  indices  for  the  two  media  and  0  and  0,  are  the  angles 
between  the  incident  and  refracted  rays  and  the  surface  normal  (see  Figure  1).  For 
incident  angles,  #,  near  grazing,  this  ratio  becomes  very  small.  As  a  consequence,  the 
apparent  angle  subtended  by  the  synthetic  aperture  becomes  very  small,  and  the  resolu¬ 
tion  in  the  direction  normal  to  the  illumination  direction  is  severely  degraded.  This  effect 
will  be  apparent  in  the  impulse  response  examples  to  follow. 


D.  OBJECTIVE  AND  LIMITATIONS  OF  ANALYSIS 

The  objective  of  the  analysis  described  here  is  to  allow  a  rapid  assessment  of  the 
impulse  response  or  spatial  resolution  that  can  be  anticipated  from  a  wide  range  of 
possible  schemes  for  subsurface  target  illumination.  The  approach  is  subject  to  several 
approximations  and  limitations  which  must  be  considered  when  interpreting  the  results 
and  generalizing  the  method  to  additional  cases: 
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Figure  1.  Effect  of  Surface  Refraction  Is  to  Compress  Angular  Spread  of  Ray  Bundles 

•  Soil  properties — We  assume  that  the  soil  is  perfectly  homogeneous.  In 
practical  systems,  soil  inhomogeneity  will  strongly  affect  the  image  quality 
obtained  [Refs.  5,  6],  Soil  losses  and  their  frequency  dependence  are 
neglected. 

•  Surface  properties — The  soil  surface  is  assumed  to  be  perfectly  planar.  We 
account  for  Snell’s  refraction  law  but  not  for  reflectance  losses  or  their 
frequency  dependence. 

•  Target-environment  interaction — Scattering  effects  within  the  target 
(resonance)  and  between  target  and  surface  are  not  considered.  Likewise, 
multiple  surface  reflections  are  not  accounted  for. 

•  Parallax  and  near-field  effects — Plane -wave  illumination  is  assumed.  We 
assume  that  the  point  at  which  illumination  rays  pass  through  the  soil  is 
(nearly)  the  same  for  all  points  in  the  synthetic  aperture.  The  effect  of  surface 
refraction  is  included  by  transforming  the  free-space  incident  angle  to  the 
refracted  angle  in  the  medium  according  to  Eq.  4. 

•  Impulse  response — Because  the  impulse  response  is  computed  within  the  soil 
medium,  we  assume  that  the  burial  depth  is  greater  than  the  depth  resolution 
of  the  system.  The  surface  does  not  appear  in  the  plots  of  impulse  response. 

•  Antenna  beamwidth — Rolloff  of  the  primary  element’s  antenna  pattern  is  not 
considered. 

•  Sidelobes  and  aperture  weighting — There  are  many  ways  in  which  target 
illumination  can  be  apodized  to  minimize  sidelobe  effects  and  account  for  the 
frequency  dependence  of  propagation  losses.  For  this  analysis,  however,  we 
consider  only  a  direct  unweighted  summation  of  the  points  in  k-space 
representing  the  target  illumination.  This  method  will  yield  a  mainlobe  width 
close  to  the  best  (minimum)  that  could  be  achieved,  but  not  necessarily  the 
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best  sidelobes.  For  this  reason  we  plot  the  impulse  response  on  a  linear  scale 
to  emphasize  mainlobe  features. 

In  most  cases  of  interest  the  results  of  this  method  will  be  good  approximations  to 
more  exact  methods,  while  the  speed  and  efficiency  afforded  allow  many  different  cases 
to  be  quickly  analyzed.  Besides  comparing  the  impulse  response  of  a  wide  range  of 
different  aperture  configurations,  we  can  evaluate  the  change  in  impulse  response  for  a 
given  synthetic  aperture  as  the  focused  pixel  is  moved  to  various  positions  in  the  radar 
field  of  regard. 
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2.  SYNTHETIC  APERTURE  ALTERNATIVES 


The  application  of  radar  to  the  systematic  detection  of  subsurface  targets, 
particularly  mines,  has  led  to  the  development  of  a  wide  range  of  alternative  systems  for 
creating  synthetic  apertures  for  target  illumination.  Early  diagnostic  systems  [Refs.  2,  3] 
scanned  a  rectangular  horizontal  aperture  directly  above  the  region  of  interest.  Later 
geometries  used  vertically  oriented  rectangular  apertures  to  examine  buried  regions 
located  in  front  of  the  radar.  Still  more  recent  systems  have  attempted  to  perform 
subsurface  imaging  with  horizontal  apertures  substantially  offset  from  the  region  of 
interest  in  a  strongly  squinted  geometry.  The  latter  two  approaches  are  particularly 
appropriate  for  forward-looking  mine  detection  since  they  could  lead  to  mine  declaration 
and  location  before  the  search  vehicle  passes  above  the  possible  target. 

Each  of  the  subsequent  four-part  figures  (Figures  2-12)  follows  a  common  plan. 
In  (A)  we  plot  first  the  real-space  region  scanned  by  the  synthetic  aperture,  with  the  point 
of  desired  image  focus  (+)  shown  for  reference.  We  then  (B)  show  the  k-space  transfer 
function  that  would  result  from  illuminating  the  focal  point  with  a  radar  bandwidth  of 
500-1,500  MHz  from  each  of  the  points  in  the  apertures.  We  then  plot  (C),  a  one¬ 
dimensional  trace  of  the  impulse  response  in  the  cross-track  direction.  This  subfigure  also 
summarizes  the  conditions  pertaining  to  the  entire  figure.  Finally,  we  show  (D),  a  two- 
dimensional  section  of  the  impulse  response  in  a  vertical  plane  through  the  focal  point. 
This  allows  the  shape  of  the  impulse  response  lobe  to  be  shown  and  also  shows  the  effect 
of  surface  refraction  in  causing  a  loss  of  resolution  in  the  direction  normal  to  that  of 
primary  illumination. 

Figure  2  shows  the  impulse  response  associated  with  a  “rail-SAR”  of  the  type 
used  by  ERIM  and  Lincoln  Laboratory  [Ref.  7]  for  the  measurement  of  both  surface  and 
subsurface  targets.  In  this  case  Figure  2A  shows  that  the  synthetic  aperture  comprises  a 
region  10  m  wide  by  5  m  high,  centered  5  m  above  the  surface.  Data  collected  from  this 
aperture  are  used  to  focus  an  image  pixel  on  the  surface  at  a  range  of  20  m  from  the 
aperture.  From  the  standpoint  of  this  pixel,  the  k-space  illumination  function  is  shown  in 
Figure  2B.  Figure  2C  is  a  one-dimensional  transverse  cross-range  trace  through  the 
amplitude  of  the  impulse  response.  Figure  2D  is  a  two-dimensional  contour  plot  of  the 
impulse  response  amplitude  in  the  vertical  plane  whose  horizontal  axis  is  the  downrange 
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direction  from  the  radar  toward  the  focused  pixel.  In  this  and  subsequent  figures,  the  five 
contour  lines  correspond  to  the  peak  of  the  impulse  response  times  (0.9,  0.7,  0.5,  0.3, 
0.1).  Figure  3  is  the  corresponding  case  in  which  the  focused  pixel  is  below  the  surface, 
so  that  refraction  has  compressed  the  aspect  angle  range  of  the  target.  Comparison  of 
Figure  2B  with  3B  shows  how  the  k-space  spectral  support  region  has  been  reoriented, 
extended  in  length  but  also  compressed  in  the  vertical  plane.  Cross-range  resolution  is 
unchanged  in  this  geometry  since  the  wave  vector  increase  in  the  dense  medium  is  offset 
by  the  angle  compression  due  to  refraction.  In  the  elevation  plane,  Figure  3D,  the  effect 
of  refraction  in  the  dense  medium  is  clear.  Comparison  with  Figure  2D  shows  that  the 
impulse  response  has  been  rotated  because  of  the  larger  angle  of  refraction,  thinned  in  the 
direction  normal  to  propagation  by  the  larger  wave  vectors,  but  extended  dramatically  in 
the  perpendicular  direction  due  to  angular  compression. 

Figures  4-7  show  a  different  synthetic  aperture  geometry,  one  typical  of  many 
ground-penetrating  radars  (GPRs)  designed  for  forward-looking  mine  detection  [Ref.  8], 
A  typical  configuration  of  these  systems  is  a  horizontal  antenna  aperture  (either  real  or 
synthetic)  extending  the  width  of  the  vehicle.  As  the  vehicle  moves  down  the  road,  radar 
sweeps  in  successive  locations  create  a  synthetic  aperture  in  a  horizontal  plane.  Because 
the  antennas  are  directed  forward,  this  amounts  to  a  strongly  squinted  synthetic  aperture 
that  can  be  focused  to  surface  and  subsurface  points  at  some  distance  in  front  of  the 
vehicle.  Figure  4  corresponds  to  a  4-m  wide  antenna  array  at  1  m  above  ground  level 
(AGL),  where  it  might  be  placed  on  a  rack  projecting  in  front  of  the  vehicle.  As  the 
vehicle  moves,  this  array  forms  a  horizontal  rectangular  (4  m  by  10  m)  synthetic  aperture 
focused  15  m  in  front  of  the  array.  Figure  4A  shows  the  real-space  aperture  and 
Figure  4B  the  corresponding  k-space  support  after  refraction.  Figure  4D  shows  that 
despite  the  substantial  increase  in  wave  vector  below  the  surface,  angle  compression  has 
essentially  negated  any  transverse  resolution  in  the  vertical  plane. 

If  true  three-dimensional  resolution  is  required,  system  designers  have  the  option 
of  increasing  the  height  of  the  antenna  array.  Figure  5  shows  the  effect  of  increasing 
aperture  height  to  6  m  AGL.  This  provides  a  slight  but  measurable  improvement  in 
vertical-plane  resolution.  Figures  6  and  7  show  the  effect  of  moving  the  processed  point 
closer  to  the  end  of  the  aperture.  This  improves  cross-range  resolution  in  both  the  vertical 
and  horizontal  planes,  although  good  three-dimensional  resolution  is  only  approached 
when  the  focused  pixel  is  adjacent  to  the  synthetic  aperture,  that  is,  when  the  standoff 
distance  is  near  zero. 
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A.  synthetic  aperture  B.  k-space  spectrum 
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Figure  2.  Rail  SAR  Impulse  Response.  Focused  point  is  in  air.  Contours  in  D.  are  at  0.9, 
0.7,  0.5,  0.3,  and  0.1  times  peak  of  impulse  response  magnitude  (volts). 
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Figure  3.  Rail  SAR  Focused  in  Dense  (n  =  2)  Medium.  Note  the  compression  of  ray  bundles 
and  loss  of  transverse  resolution  in  the  elevation  plane. 

The  conclusion  from  Figures  4-7  is  that  it  will  be  very  difficult  for  forward- 
looking  GPR  systems  to  achieve  both  significant  subsurface  vertical  resolution  and  good 
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standoff  range.  For  mine  detection,  this  may  be  acceptable,  provided  significant  sub¬ 
surface  resolution  is  not  needed  to  reject  clutter. 


C.  crossrange  ipr 


B.  k-space  spectrum 


0  -4@ 


ky/(rad/m)  -  -  kx/(rad/m) 

D.  downrange-vertical  ipr 


Figure  4.  Forming  a  Synthetic  Aperture  by  Forward  Motion  of  a  Transverse  Linear  Array. 
Near-grazing  angle  leads  to  total  compression  of  angular  spread  and  loss  of  resolution. 
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Figure  5.  Raising  the  Antenna  Array  Provides  a  Larger  Grazing  Angle.  Transverse  vertical 

plane  resolution  is  slightly  improved. 
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A.  synthetic  aperture 


B.  k-space  spectrum 
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Figure  6.  Both  Raising  the  Antenna  and  Forming  an  Image  Closer  to  the  Synthetic 
Aperture  Shows  Significant  Improvement  in  Vertical  Plane  Resolution 
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Figure  7.  The  Extreme  Case  Forms  the  Image  at  the  Edge  of  the  Synthetic  Aperture.  This 
provides  the  best  resolution  but  at  minimal  standoff  range. 
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A  different  synthetic  aperture  configuration,  which  has  been  widely  used  in 
diagnostic  GPR,  is  a  rectangle  placed  directly  over  the  region  of  interest.  This  maximizes 
the  symmetry  of  the  k-space  transfer  function  support  and  also  provides  what  is,  in 
principle,  the  largest  support  affordable  in  a  given  frequency  band.  In  Figure  8  we 
analyze  the  response  of  a  2-m  square  aperture  centered  above  the  pixel  of  interest.  In  this 
case  the  focal  point  is  assumed  to  be  in  free-space  ( n  =  1).  Note  from  Figure  8D  that  the 
downrange  resolution  is  about  twice  as  fine  as  the  vertical  resolution.  This  is  consistent 
with  Figure  8B,  which  shows  that  the  downrange  and  cross-range  support  of  the  transfer 
function  is  about  twice  its  vertical  extent.  The  geometry  is  the  same  in  Figure  9,  except 
that  the  focal  point  is  in  a  soil  medium  ( n  =  2).  Note  that  while  the  vertical  resolution  is 
twice  as  fine  as  in  Figure  8D,  transverse  (horizontal)  resolution  is  the  same  as  in  the 
previous  free-space  case.  This  result  is  a  consequence  of  Eqs.  3  and  4.  Immersing  the 
focal  point  in  the  medium  of  index,  n,  compresses  sin  (A0/2)  and  X  by  the  same  factor  of 
n,  leaving  transverse  resolution  unchanged.  Figure  10  is  similar  to  Figure  9,  except  that 
the  focused  pixel  is  offset  0.7  m  from  the  center  of  the  aperture.  The  effect  of  asymmetric 
illumination  (Figure  10D)  is  to  change  the  sidelobe  structure  without  significantly 
affecting  mainlobe  resolution.  This  suggests  that  a  single  synthetic  aperture  position  can 
provide  a  large  field  of  regard  for  imaging,  without  unacceptable  loss  of  image  quality. 

Because  good  resolution  in  three  dimensions  arises  not  only  from  wideband 
illumination  but  also  from  the  widest  possible  range  of  diversity  in  illumination  direction, 
it  is  worthwhile  to  consider  the  most  efficient  aperture  shapes  from  the  standpoint  of 
maximizing  resolution  from  the  smallest  number  of  measurement  points.  One  possible 
choice  for  an  efficient  synthetic  aperture  is  a  circular  trajectory  [Ref.  9],  as  shown  in 
Figure  11  A.  This  corresponds  to  moving  a  radar  20  m  AGL  in  a  circle  of  20-m  radius 
about  the  focused  pixel.  Figure  11B  shows  the  corresponding  k-space  support  for 
subsurface  focusing.  Note  the  overall  similarity  in  the  support  volume  between  this  case 
and  that  of  the  filled  square  aperture  described  in  Figure  9B.  Surface  refraction  tends  not 
only  to  move  the  incidence  angle  of  the  illumination  rays  closer  to  the  vertical,  but  also  to 
compress  their  range  of  variation.  As  a  consequence,  there  is  less  difference  in  k-space 
support  between  the  filled  square  and  hollow  circular  apertures  than  would  at  first  be 
expected.  This  is  further  confirmed  by  comparison  of  the  impulse  responses  of  the  two 
cases  (Figures  9D  and  1  ID).  Despite  the  apparent  differences  in  illumination,  the 
mainlobe  shapes  of  the  impulse  responses  are  very  similar. 
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A.  synthetic  aperture 


B.  k-space  spectrum 


ky/(rad/m)  -20  -20  kx/(rad/m) 


C.  crossrange  ipr  D.  downrange-vertical  ipr 


Figure  8.  Typical  Geometry  for  Tomographic  or  Holographic  Reconstruction. 
This  pixel  is  imaged  in  air  (n  =1). 


A.  synthetic  aperture 
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C.  crossrange  ipr 


D.  downrange-vertical  ipr 


Figure  9.  Same  Geometry  as  Figure  8,  but  Pixel  Is  Focused  in  Dense  Medium.  Vertical 
resolution  is  improved  but  transverse  is  unchanged. 
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A.  synthetic  aperture 


B.  k-space  spectrum 


C.  crossrange  ipr 


D.  downrange-vertical  ipr 


horizontal  distance/meters  downrange/m 


Figure  10.  Focusing  an  Off-Center  Point  Leads  to  Some  Distortion  of  Impulse  Response, 
but  Volume  of  Resolution  Cell  Is  Negligibly  Affected 


A.  synthetic  aperture 


B.  k-space  spectrum 
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Figure  11.  The  Empty-Circle  Aperture  Provides  Similar  Illumination 
Angle  Diversity  to  Figure  9.  Note  similarity  of  resolution  cell. 
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The  final  point  to  be  made  is  that  the  circular  aperture  is  similarly  robust  with 
respect  to  off-center  focusing.  In  Figure  12  we  consider  the  case  of  a  focused  pixel  10  m 
from  the  aperture  center.  As  we  saw  for  the  case  of  the  filled  square  aperture  (Figures  9 
and  10),  although  the  shape  of  the  impulse  response  mainlobe  and  sidelobe  structure  is 
changed,  the  basic  resolution  remains  the  same  as  for  the  centered  case. 


A.  synthetic  aperture 


B.  k-space  spectrum 


C.  crossrange  ipr 


D.  downrange-vertical  ipr 


horizontal  distance/meters  downrange/m 


Figure  12.  Focusing  an  Off-Center  Point  Changes  Sidelobe  Structure  with  Negligible 

Impact  to  Mainlobe  Resolution 


13 


3.  CONCLUSIONS 


Ground-penetrating  radar  is  an  important  technology  for  subsurface  target 
detection,  with  particular  application  to  both  mine  and  UXO  detection.  Although  it  has 
been  in  development  for  many  decades,  it  is  clear  that  much  work  remains  to  be  done. 
Specific  designs  for  imaging  subsurface  radar  have  been  proposed  by  many  organiza¬ 
tions,  and  several  have  undergone  significant  system  development.  Before  making  major 
design  commitments,  however,  it  is  desirable  to  be  able  to  estimate  the  anticipated 
resolution  of  a  given  system  concept  based  on  a  first-order  map  of  the  k-space 
illumination  support  its  geometry  will  afford.  The  techniques  described  in  this  paper 
allow  a  system  designer  to  explore  the  effect  of  variations  in  system  mechanical  design 
on  radar  resolution.  In  addition,  we  are  able  to  estimate  the  field  of  regard  which  may  be 
focused  from  a  single  aperture  without  excessive  degradation  of  image  quality.  Finally,  it 
is  clear  that  the  angular  spread  of  ray  bundles  is  severely  compressed  by  surface 
refraction,  particularly  for  near-grazing  incidence,  with  consequent  impact  on  subsurface 
resolution. 
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